Controlling spatiotemporal chaos 
in oscillatory react ion- diffusion systems 
by time-delay autosynchronization 



C. Beta, A. S. Mikhailov*, 

Ahteilung Physikalische Chemie, 
Pritz-Haber-Institut der Max-Plack-Gesellschaft, 
Faradayweg 4^6, 14-195 Berlin, Germany 



Abstract 

Diffusion-induced turbulence in spatially extended oscillatory media near a super- 
critical Hopf bifurcation can be controlled by applying global time-delay autosyn- 
chronization. We consider the complex Ginzburg-Landau equation in the Benjamin- 
Feir unstable regime and analytically investigate the stability of uniform oscillations 
depending on the feedback parameters. We show that a noninvasive stabilization 
of uniform oscillations is not possible in this type of systems. The synchronization 
diagram in the plane spanned by the feedback parameters is derived. Numerical 
simulations confirm the analytical results and give additional information on the 
spatiotemporal dynamics of the system close to complete synchronization. 
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1 Introduction 



Pattern formation in spatially extended systems far from thermal equilibrium 
has attracted much attention in recent years and became a field of active 
research [1,2]. A large class of these systems is described by equations of the 
reaction-diffusion type. In such systems, a broad variety of complex spatiotem- 
poral patterns has been observed, ranging from travelling pulses and target 
patterns to rotating spiral waves and spatiotemporal chaos. Important exam- 
ples are chemical reaction-diffusion systems such as the Belousov-Zhabotinsky 
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reaction or catalytic surface chemical reactions [3,4,5], as well as semiconduc- 
tor and biological systems [6,7]. 

Many systems show oscillatory dynamics of their individual elements. Close to 
the onset of oscillations, the equations describing any oscillatory system can 
be reduced to a simple universal equation for the complex oscillation ampli- 
tude, the complex Ginzhurg- Landau equation (CGLE) [8,9]. For a distributed 
system with diffusive coupling, the CGLE corresponds to the normal form of a 
supercritical Hopf bifurcation. Although the CGLE is strictly vahd only in the 
immediate vicinity of the oscillation onset, it has been found in many cases 
that its qualitative predictions also hold within a wider range near the bifurca- 
tion point. The CGLE thus provides a general framework for studying common 
features of the dynamical behaviour in oscillatory reaction-diffusion systems. 
If oscillations are synchronized, uniform oscillations can be estabhshed in the 
medium or regular spatiotemporal patterns can emerge. 

If the Benjamin- Feir criterion is met, synchronous oscillations become unsta- 
ble and a turbulent state develops where two types of chaotic behaviour can be 
distinguished. In the state of phase turbulence, the local oscillation phase ex- 
hibits weak irregular fluctuations while the real amplitude remains saturated 
showing only small variations [8]. Amplitude or defect turbulence is charac- 
terized by strong fluctuations of both phase and amplitude which are due to 
the presence of defects, the locations with a vanishing amplitude [10]. 

Over the past decade, problems of chaos control became to play a central role 
in the stTidics of chaotic dynamics [11]. Inspired by the pioneering work of 
Ott, Grebogi, and Yorke (OGY) [12], control of chaos has been studied in the 
context of many different dynamical systems [13,14,15] However, the OGY 
method is designed to control chaos in low-dimensional dynamical systems. It 
requires continuous extensive analysis of the system dynamics which is virtu- 
ally impossible to carry out in the case of fully developed high-dimensional 
turbulence in spatially extended systems. A different empirical control method 
based on implementing a delayed feedback loop was proposed by Pyragas and 
is known as time-delay auto synchronization (TDAS) [16,17]. In this method, a 
feedback signal F is apphed to the system that is proportional to the difference 
between the delayed value of a given system variable u and its instantaneous 
value, F ~ u{t — r) — u{t). This approach is easily extended to spatially dis- 
tributed experimental systems that often do not allow individual addressing 
of their local elements. Then, the global feedback signal is generated from 
the integral value of u over all system elements. In recent studies, the TDAS 
method has been used in a number of different dynamical systems [18,19,20] 
and was modifled and improved considerably [21,22,23]. For the CGLE, the 
action of local TDAS has been previously considered [24,25]. 

The present work was motivated by our experimental study of catalytic CO 
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oxidation under TDAS [26]. The catalytic CO oxidation on a platinum (110) 
single crystal surface constitutes the most thoroughly studied example among 
a class of simple heterogeneous catalytic reactions for which reaction-induced 
spatiotemporal pattern formation has been reported [5]. It displays a partic- 
ularly rich variety of spatiotemporal concentration patterns [27] that can be 
modeled numerically using a simple model of three coupled partial differential 
equations [28,29]. The effect of global delayed feedback on pattern formation 
in catalytic CO oxidation was studied previously using a different time-delay 
feedback scheme [30,31,32]. The same feedback scheme was also investigated 
in detail in the context of the complex Ginzburg-Landau equation [33,34]. 
Our work on suppression of chemical turbulence in catalytic CO oxidation 
using the TDAS scheme [26] was focused on investigating the question of in- 
vasiveness. We found in the experiment that, although we could reduce the 
magnitude of the feedback signal considerably, we were not able to reach the 
ideal limit of a vanishing feedback signal in the state of control. In the present 
work, we use the complex Ginzburg-Landau equation as a general model for 
spatially extended oscillatory systems to study common aspects of control of 
diffusion-induced chemical turbulence by TDAS. 

The paper is organized as follows. In the next section, the physical problem 
and the corresponding equations will be introduced. In Section 3, we analyze 
the stability of uniform oscillations in the presence of TDAS and derive a 
synchronization diagram. Numerical simulations of the complex Ginzburg- 
Landau equation with TDAS are described in Section 4. The paper ends with 
conclusions and a discussion of the obtained results. 



2 Formulation of the problem 

In absence of feedback, the dynamics of a chemical reaction-diffusion system 
can be generally described by a set of coupled partial differential equations 



where u = {ui,U2, ■■■,uiy) represents concentrations of reacting species and D 
is their diffusion matrix. The set of nonlinear functions f (u, p) of the concen- 
trations u accounts for the reaction part of the dynamics and depends on a 
number of parameters p = (j>i,P2, ■■■,Pm), such as rate constants and external 
conditions. 

Suppose that the dynamics of system (1) is such that a supercritical Hopf 
bifurcation occurs when crossing a certain threshold p = po in the parameter 
space. Close to the onset of oscillations, system (1) can be transformed to a 



u = f(u,p) + DV2u, 
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simple description of the dynamics in terms of a complex oscillation amplitude 
rj{x,t). This transformation is based on retaining only the leading critical 
modes that govern the dynamics close to the bifurcation point. It yields as 
equation of motion for the amplitude ri{x,t), the complex Ginzburg-Landau 
equation [8], 

?7 = (1 - iuj)r] - (1 + i/3) \r]\^r] + {l + is) V^t? , (2) 

where the parameters u, (3, and e denote the linear oscillation frequency, the 
nonlinear frequency shift, and the linear dispersion coefficient, respectively. 
When the condition 1 + e/3 < of the Benjamin-Feir instability is satisfied, 
uniform oscillations are unstable and turbulence spontaneously develops in the 
system. To control the turbulence, global feedback can be introduced. In the 
presence of global time-delayed feedback, the considered system is described 
by the equation 

7) = (1 - iu;)r] - (1 + i/3) 1771^77 + (1 + i£) V'?7 + F(t) , (3) 

where the feedback term is given by 

F{t)=i,e'^rj{t-T)-rj{t)), (4) 

and 

1 f 

V{t) = r]{x,t)dx. (5) 



Here, /i is the feedback intensity factor, r is the delay time, and x is a phase 
shift in the application of the control force. 

Previously, we showed experimentally [26] that diffusion-induced turbulence 
could be efficiently suppressed in an oscillatory chemical reaction-diffusion 
system using TDAS. It was found that the invasiveness of the control scheme 
could be significantly reduced for an optimized choice of feedback parameters. 
Figure 1 presents the period of uniform oscillations in the state of control and 
the feedback magnitude as a function of the delay time in the experiments with 
CO oxidation on Pt(llO). Although the feedback magnitude was significantly 
reduced for a delay close to r = 6 s, the optimal case of a vanishingly small 
feedback signal could not be reached. Near the point, for which the period T of 
stabilized uniform oscillations becomes equal to the delay time r, an instability 
was found. As seen in Fig. 1, the system avoids a state corresponding to the 
intersection point with the line for which T = r. The stability of this state, for 
which the feedback vanishes, has been earlier investigated [26] for the uniform 
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system in terms of a phase dynamics equation of a single oscillator. In this 
article, we extend the analysis to the theoretical study of a spatially extended 
system. 

Since our present work was motivated by the experimental results in [26], we 
restrict our investigation to the case of a globally applied control signal gen- 
erated from the averaged complex amplitude. This case is relevant in many 
experimental situations where a controlled manipulation of the individual sys- 
tem elements is not possible. Note, however, that a different behaviour can be 
expected for a space dependent application of the control scheme. 



3 Linear stability analysis 

We choose a superposition of a homogeneous mode H with a small spatially 
inhomogeneous perturbation of arbitrary non-zero wave number k as an ansatz 
for the complex oscillation amplitude, 

r){x, t) = H{t) + A+{t)e'^'^ + A_{t)e-"''^ . (6) 

Substituting expression (6) into Eq. (3), we can separate, for small A±, homo- 
geneous contributions from the spatially inhomogeneous terms. In this way, 
we find an equation for the homogeneous mode which is decoupled from the 
nonuniform contributions, 

if = (1 - ioj)H - (1 + i(3)\H\'^H + ^jie'''{H{t - r) - H(t)) , (7) 

and a pair of coupled equations for the wave amplitudes A±, 

i+= (1 - iuj)A+ - {l + \e)K^A+ - 2{l+iP)\H\^A+ - {l + iP)H^A*_, (8) 
A*_ = (1 + iu)A*_ - (1 - ie)K^A*_ - 2(1- il3)\H\''A*_ - (1 - il3)H*''A+.{Q) 

Let us first discuss the dynamics of the uniform mode described by Eq. (7). 
Substituting H = poe"^"* into Eq. (7) we obtain the amplitude of uniform 
oscillations in the presence of TDAS, 

Po = ^/l + /^(cos(x + ^t) - cosx) , (10) 

and derive the equation 

— uj -\- [3 + nP{cos{x + ^t) — cosx) ~ /^(sin(x + ^t) — sinx) (H) 
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determining the frequency Q of uniform oscillations. Note that without the 
feedback (/x = 0) the system shows uniform oscillations with the frequency 

no = '^ + P- 

Figure 2 shows the results of numerical integration of Eq. (7). The oscillation 

period T = 27r/fl is shown as a function of the delay time r for weak and strong 
feedbacks. We see that the state with T = r is stable for a weak feedback 
and becomes unstable for a high feedback intensity. The strong feedback case 
[Fig. 2(b)] shows a good qualitative agreement with the experimental data in 
Fig. 1. 

The stability of the state with T = r can be understood in terms of the 
bifurcation diagram presented in Fig. 3. This diagram is constructed by solving 
Eq. (11) for a delay time r equal to the period Tq = 2n/{uj+(3) oi oscillations in 
the nonperturbed uniform system. Besides a solution with a vanishing feedback 
term, for which = r^oi other solutions with a non-zero feedback and Q ^ 
are yielded by Eq. (11). Linear stability analysis shows that the solution with 
r2 = and a vanishing feedback is stable for small /j, and becomes unstable 
beyond some critical feedback intensity Hq. These results are similar to those 
obtained by using the phase dynamics approximation for a single oscillator in 
the presence of TDAS [26]. 

We now turn to the stability analysis of uniform oscillations with respect to 
spatially inhomogeneous perturbations. The solution of the linear equations 
(8) and (9) can be sought as = A° e-*"*e^* and A*_ = Al°e^"*e^*. Substi- 
tuting this into (8) and (9), we obtain the following expression 

Ai,2 = 1-k'-2pI± ^(1 + (3^)pt -{n-u-SK^- 2/3piy , (12) 



where po is given by Eq. (10) and Q is a solution of Eq. (11). Uniform oscil- 
lations are stable with respect to the growth of spatially nonuniform modes 
if Re Ai,2 < for all wavenumbers k. The instability boundary is therefore 
determined by the conditions 



ReA(//c, «c) = 0, (13) 
^ReA(//e,«e) = 0. (14) 

Here, Eq. (14) accounts for the fact that the wave number of the first un- 
stable mode corresponds to a maximum in A plotted as a function of n. The 
equations (13) and (14) can be solved numerically, taking into account Eq. 
(11) yielding the dependence of the frequency Q of the uniform mode on the 
feedback intensity p and the delay time r. 
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We are particularly interested in the stability of uniform oscillations with 
^2 = ^0 and T — Tq. Inserting this into the general expression (12) for A, we 
find 



Al,2 = -K' 



(15) 



which is independent of /x. Therefore, this state is always unstable when the 
Benjamin- Feir condition l-|-£/3 < is fulfilled. All inhomogeneous modes with 
a wavenumber less than 



are growing, no matter how jj, is chosen. Thus, for 1 + £/5 < the solution 
with — cu + (3 will be always unstable so that a noninvasive stabilization of 
uniform oscillations with TDAS is not possible in this type of system. 

For the other solutions with fl ^ Qq at t = Tq, presented in Fig. 3, the feedback 
signal is nonzero. In this case, we can determine fj,c from the general expression 
(12) for A by numerically solving equations (13) and (14). Figure 4 shows the 
critical feedback intensity as a function of the dispersion parameter e for 
T = Tq and the other parameters chosen as in Fig. 3. In the Benjamin-Feir 
unstable regime, /ic lies above the bifurcation point Hq at which the solution 
with ft — flQ becomes unstable, cf. Fig. 3. As the system approaches the 
Benjamin-Feir line with decreasing e, the critical feedback intensity necessary 
to stabilize uniform oscillations decreases and finally converges towards /iq. 

From the conditions (13) and (14), we can also determine the critical feedback 
intensity /^c as a function of the delay time r. The resulting synchronization 
diagram in the plane spanned by the feedback parameters r and /i is displayed 
in Fig. 5(a). The curve for the critical feedback intensity Hc divides the plane 
into a shaded region, where uniform oscillations are linearly stable with respect 
to small perturbations of arbitrary wavenumber, and a region where uniform 
osciallations are unstable. The characteristic feature of this diagram is the 
repeated appearance of cusps. They are observed whenever r becomes equal 
to an integer multiple of the period of the unperturbed uniform system, r = 
k27r/{uj + f3),k = 1,2,3,... . Figure 6 shows an extension of the top part 
of Fig. 5 towards large delay times. With increasing r, the cusps get less 
pronounced and the boundary converges to a flat line at /i ~ 0.16. 

According to this analytically derived synchronization diagram, stability of 
uniform oscillations can also be maintained by applying global feedbacks with 
very large delay times. Moreover, the critical feedback strength, needed to 
maintain synchronization, does not depend on the delay in the limit r — > oo. 




(16) 



7 



To understand this result, we note that the feedback signal F{t) for r 
is given by 



oo 



F{t) = fie 



^X 



Pod 



-i{nt+(t>o) 



(17) 



where 00 is a constant phase shift. The first term here corresponds to the state 
f]{t — t) at t — > oo which is essentially the initial state of the system. When 
destabilization of initally uniform oscillations is considered to determine the 
stability boundary, this initial state represents uniform oscillations with fre- 
quency Q and amplitude po- Substituting this expression for F(t) into Eq. (3), 
we see that then a situation with external periodic forcing is effectively real- 
ized. The critical value pic corresponds in this case to the minimum forcing 
intensity needed to maintain uniform oscillations in the system. 

In Fig. 5(b,c), the critical wavcnumber and the frequency arc shown, respec- 
tively, as functions of the delay time along the lower part ABC of the stability 
boundary. The two curves are of similar shape: they display a decrease for in- 
creasing delay time interrupted by discontinuous jumps. These discontinuities 
occur at the locations where the cusps are found in the stability boundary in 
Fig. 5(a). 

In Fig. 7, we show the real part of A as a function of k, at three different 
points on the stability curve in the (p, r) plane. The three cases correspond 
to what has been found earlier for a different global delayed feedback scheme 
in the CGLE [33]. On the branch AB (and similarly also on the branches BC 
and to the right of C), the first unstable modes occur with a wavenumber 
Kq where Im(A) = 0, see Fig. 7(b,c). Thus, if we cross this branch of 
the stability boundary by reducing the feedback intensity below the critical 
value yUc, uniform oscillations will become unstable and standing waves with 
wavenumber kq will emerge (cf. Sec. 4). A different situation is encountered on 
the branch reaching from A upwards, see Fig. 7(a). Here, the instability will 
occur by periodic spatiotemporal modulations of uniform oscillations, since 
Im(A) 7^ and the most unstable modes will have wavenumbers close to 
Kq = 0. 



4 Numerical simulations 



All simulations were carried out for a one-dimensional system of length L — 
128 on an equidistant grid with Ax = 0.32 (corresponding to a total num- 
ber of 400 grid points). We use a second-order finite-difference scheme for the 
discretization of the Laplacian operator and impose periodic boundary con- 
ditions. An explicit Euler scheme with a fixed time step of At — 0.001 was 
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employed for integration. As initial condition, a uniform state sTipcrposed with 
a small spatially inhomogeneous perturbation was chosen. The set of param- 
eters is as in the previous section: e = 2, /3 = —1.4, lu = 27t — (3 ^ 7.68, and 
X = vr/2. The choice of the feedback parameters // and r is different for the 
different simulations and specified below. 

In the first series of simulations, we systematically scan the parameters fi and r 
in steps of A/j, = At = 0.1 between and 2.5, respectively, to verify the shape 
of the stability domain in the {/i, r) plane. When starting from uniform initial 
conditions, the stability diagram displayed in Fig. 5(a) is nicely reproduced 
after transients. Asymptotic states are uniform inside the domain of stability 
and nonuniform outside this region. These nonuniform states are, however, of 
different types. Far from the region of stability, we asymptotically reach a fully 
developed state of defect-mediated turbulence. As the stabihty boundary is 
approached, phase turbulence and standing wave patterns are observed close 
to the branches AB, BC, and to the right of C. 

Figure 8 shows the results of a more detailed scan of n and r in the vicinity 
of the branch AB of the synchronization diagram in Fig. 5(a). Again, we start 
from uniform initial conditions. The feedback parameters n and r are changed 
in steps of A/x = At = 0.025 between r = 0.25 ... 0.975 and = ... 0.7. Simu- 
lations that, after transients, resulted in a regular nonuniform spatiotemporal 
state are marked by bold dots at the corresponding (/i, r) coordinates (sim- 
ulations leading to a uniformly oscillating or turbulent asymptotic states are 
not shown). Obviously, the parameter region, where regular spatiotemporal 
patterns occur, constitutes a slightly asymmetric prolongation of the tongue- 
shaped stability domain for uniform oscillations towards smaller feedback in- 
tensities. From the result of the first coarse parameter scan we conjecture, that 
the regions where spatiotemporal patterns occur look qualitatively similar at 
the other tongue-shaped branches of the stabihty curve. 

When starting from turbulent initial conditions, the stability boundary for 
uniform oscillations is moved towards larger feedback intensities in the vicin- 
ity of the cusps. The results of numerical simulations are summarized in Fig. 9. 
Here, open circles indicate a turbulent asymptotic state and bold circles again 
denote the appearance of a regular nonuniform wave pattern. Inside the area 
where no symbols are shown, simulations converge towards uniform oscilla- 
tions. For comparison, the analytically derived synchronization diagram from 
Fig. 5(a) is also displayed. 

To get a better idea of how the spatiotemporal dynamics of the system changes 
in the parameter range between turbulence and uniform oscillations, we show 
in Fig. 10 a series of spacetime plots of asymptotic dynamical states reached for 
different feedback intensities at a fixed delay time of r = 0.5. The amplitude |?7| 
is plotted in a grey scale color coding and the feedback intensity /i is increased 
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in Fig. 10 from (a) to (c). In absence of feedback (Fig. 10(a), /i = 0) and for 
small feedback intensities (Fig. 10(b), fi = 0.05), an irregular state of defect- 
mediated turbulence is observed. However, the number of defects is smaller in 
the presence of a weak feedback and the time evolution shows intervalls, where 
almost no defects are seen. If the feedback intensity is increased (Fig. 10(c), 
H = 0.07), defects are no longer observed and the system displays a disordered 
state of phase turbulence. For even stronger feedback, breathing (Fig. 10(d), 
/J, = 0.1) and stationary standing wave patterns (Fig. 10(e), = 0.15) can be 
observed. For still larger feedbacks, /i > 0.2 (not shown), uniform oscillations 
take place. 

To study the behaviour of the system for very large delays, a series of simula- 
tions with the delays r varying from 0.5 to 29.5 in steps of At = 1 has further 
been performed. The feedback intensity ^ was varied for each choice of r from 
0.1 to 0.3 in the steps of Afj, = 0.1. The initial conditions in these simulations 
represented the state of amplitude turbulence reached by the system if /i = 0. 
These numerical investigations have revealed that synchronization is possi- 
ble for all chosen delays, when the feedback intensity exceeds a critical level. 
In another series of experiments, the feedback intensity was fixed at /i — 1 
and the delay time r was increased from to 50 in steps of 0.1. Here, the 
transient time was determined for each of the simulations by looking at the 
convergence of the statistical variance of the amplitude p. It was found that 
for small delays the transient time increases proportionally to the delay r. For 
large delays, the transient time undergoes saturation at a level independent 
of T, but depending on the intensity of the applied feedback. 

For very long delays r, the component f]{t — t) in the feedback signal F{t) 
at time t corresponds to the initial state of amplitude turbulence. Thus, the 
global delayed feedback scheme in the considered limit effectively represents 
global forcing of the system with an external chaotic signal corresponding to 
the initial turbulent state in absence of feedback. To verify this conjecture, 
special numerical simulations have been performed. In these simulations, the 
feedback signal was generated by replacing f}{t — t) with the average complex 
oscillation amplitude of the same system without the feedback. We have found 
that, by applying such chaotic external forcing, synchronization of uniform 
oscillations can also be achieved. 



5 Discussion 

Motivated by our experiments [26] on control of chemical turbulence in the 
CO oxidation reaction on Pt(llO), we have performed a detailed study of the 
effects of time-delay autosynchronization on uniform oscillations in a general 
model described by the complex Ginzburg-Landau equation. We have seen 
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that, like in the experiments, the noninvasive stabihzation of uniform osciUa- 
tions by this method is not possible, though the required magnitude of the 
feedback signal can be significantly reduced by using an optimal delay time. 
The obtained synchronization diagram exhibits a series of cusps at the values 
of delay times equal to integer multiples of the oscillation period of the unper- 
turbed uniform system. This feature seems to be common for various oscilla- 
tory systems with delayed coupling and has also been found in the case of the 
Kuramoto model of phase oscillators with a delayed global coupling [35] . Near 
the synchronization boundary, the formation of standing waves and a state of 
modified intermittent turbulence were observed. These effects are similar to 
those previously discussed for a different feedback scheme [32,33,34]. In con- 
trast to our previous studies, the action of feedbacks with long delay times 
has also been investigated. We have shown that such feedbacks are also capa- 
ble of synchronizing oscillations and this effect can essentially be explained as 
suppression of turbulence by external noisy forcing. 
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6 Figures 




Fig. 1. Experimental results on the application of TDAS to the catalytic CO oxida- 
tion on Pt(llO). Period T of homogeneous oscillations (black squares) and feedback 
magnitude M (open circles) in dependence on the delay time r. Reproduced from 
Ref. [26]. 




0.6 ^ ' ' ' ' ' ' ' 0.6 ^ ' ' ' ' ' ' ' 

0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 

delay time x delay time x 



Fig. 2. Results of numerical integration of Eq.(7). The oscillation period is plotted 

as a function of the delay time for (a) low (// = 0.5) and (b) high (/i = 1.0) feedback 
intensities. The parameters are /? = —1.4, u) = 2tt — {3 k, 7.68, and x = 
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Fig. 3. Bifurcation diagram for t = Tq = I. The parameters are as in Fig. 2. 
Dotted lines denote unstable branches. The uniform solution with 0, = = 2Tr 
and a vanishing feedback signal undergoes a transcritical bifurcation at = //q and 
becomes unstable. 
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Fig. 4. Critical feedback intensity fXc for r = Tq = 1 as a function of the dispersion 
coefficient e with (3 = —1.4. The line jj, = jiQ denotes the feedback intensity for 
which the transcritical bifurcation occurs in the uniform system (cf. Fig. 3). The 
other parameters are as in Fig. 2. 
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Fig. 5. (a) Synchronization diagram. Uniform oscillations are stable inside the 
shaded region. (b,c) The dependences of the critical wavenumber Kc and the crit- 
ical frequency Qc on the delay time r. The parameters are e = 2, (3 = —1.4, 
uj = 2Tr-P^ 7.68, and x = V^- 





Fig. 7. Growth rate Re A of spatially nonuniform modes as a function of wavenum- 
ber K at three different points on the stability curve displayed in Fig. 5(a). The 
parameters are r = 0.29, = 1.29 (a), r = 0.36, /j, = 0.4 (b), and r = 0.7, n = 0.16 
(c); other parameters as in Fig. 5. 
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Fig. 8. Regular spatiotemporal patterns at the border of complete synchronization 
for uniform initial conditions. Section of the (p, r) plane with the analytical result 
for the AB branch of the stability curve (cf. Fig 5) plotted as a solid line. Bold dots 
mark the parameters where regular spatiotemporal patters are observed. 
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Fig. 9. Summary of numerical simulations starting from turbulent initial conditions. 
Asymptotic states are shown in the {h,t) plane. Open circles denote amplitude 
turbulence and bold circles indicate regular nonuniform wave patterns. Symbols arc 
omitted in the areas where computations converge to uniform oscillations. The solid 
line shows the analytically derived stability boundary for uniform initial conditions. 
The parameters are as in Fig. 5. 
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Fig. 10. Spatiotcmporal patterns in the parameter region between turbulence and 
complete synchronization. The amplitude is displayed in a grey scale coding, 
where black (white) denotes low (high) values of the real amplitude. The delay time 
is kept constant, r = 0.5 s, and the feedback intensity increases from top to bottom, 
(a) /X = 0, (b) /X = 0.05, (c) n = 0.07, (d) n = 0.1, and (e) n = 0.15. 
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